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Abstract 

Since gene expression approaches constitute a starting point for investigating plant-patliogen systems, we performed a 
transcriptional analysis to identify a set of genes of interest in tomato plants infected with F. oxysporum f. sp. lycopersici {Fol) 
and Tomato IVlosaic Virus {ToMV). Differentially expressed tomato genes upon inoculation with Fol and ToMVwere identified 
at two days post-inoculation. A large overlap was found in differentially expressed genes throughout the two incompatible 
interactions. However, Gene Ontology enrichment analysis evidenced specific categories in both interactions. Response to 
ToMV seems more multifaceted, since more than 70 specific categories were enriched versus the 30 detected in Fol 
interaction. In particular, the virus stimulated the production of an invertase enzyme that is able to redirect the flux of 
carbohydrates, whereas Fol induced a homeostatic response to prevent the fungus from killing cells. Genomic mapping of 
transcripts suggested that specific genomic regions are involved in resistance response to pathogen. Coordinated 
machinery could play an important role in prompting the response, since 60% of pathogen receptor genes (NB-ARC-LRR, 
RLP, RLK) were differentially regulated during both interactions. Assessment of genomic gene expression patterns could 
help in building up models of mediated resistance responses. 
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introduction 

Disease resistance in plants depends on the ability of the host to 
recognize pathogens and initiate defense mechanisms that limit 
infection. A basal type of immunity in plants is conferred by the 
recognition of conserved microorganism-associated molecular 
patterns (MAMPs) by specific pattern-recognition receptors 
(PRRs) that protect hosts against non-specialized pathogens. 
Plants are also capable of establishing immune responses by using 
pathogen receptors (mainly Nucleotide binding- ARC /leucine-rich 
repeat proteins NB-ARC-LRR but also Receptor Like Proteins 
RLP; Receptor-Like Kinase RLK) known as resistance (R) 
proteins able to recognize the presence of pathogen effector 
molecules and to activate effector-triggered immunity (ETI). This 
response is typically associated with programmed cell death of the 
infected cells and the production of antimicrobial molecules [1]. 
Modulation of hormone pathways is required to restrict pathogen 
invasion, re-allocate resources, control cell death and modify plant 
architecture [2]. In addition, systemic-acquired resistance, which 
immunizes against subsequent infections, could also occur [1]. 

Tomato [Solanum Ijicopersicum) has served as an important model 
system for studying the genetics and molecular basis of resistance 
mechanisms in plants. The breadth of pathogen classes affecting 
tomato underscores the importance of tomato pathosystems as 
amenable models for studying the plant immune system [3] . Plants 



build up appropriate defense responses without draining energy 
resources to unsustainable levels through the cross-talk and fine- 
tuning of different defense pathways [4]. Identification of host 
genes involved in resistance responses could be helpful both for 
understanding plant defence mechanisms against pathogens and 
building up a biological model of plant-pathogen interaction. 

Fusarium oxysporum £ sp. lycopersici {Fol) is one of the main diseases 
affecting tomato, causing vascular wilt disease by colonizing the 
xylem vessels of roots and stems. Three physiological races, termed 
1 , 2 and 3, have been identified. The FOL A VR genes {A VRl, A VR2 
and AVR3) are carried in different combinations in different FOL 
races. AVRl ( = SIX4) is unique to race 1 whereas AVR2 ( = SLX3) is 
found in races 1 and 2. AVR3 ( = SIXl), which exists in all races /, 
L2 and L3 are known confer resistance to FOL in tomato. Race 1- 
resistant cultivars [L i2 i3), races 1 and 2-resistant cultivars {L L2 i3), 
and races 1, 2 and 3-resistant cultivars (/ L2 L3) have been 
developed [5]. To date only 1-2 gene, a coUed-coH/nucleotide 
binding- ARC /leucine-rich repeat (CC-NB-ARC-LRR) protein, 
that confers resistance to race 2 of Fusarium oxysporum f. sp. lycopersici 
has been cloned [6;7]. Responses mainly involves the callose 
deposition, the accumulation of phenolics and the formation of 
tyloses (outgrowths of xylem contact cells) and gels in the infected 
vessels [8]. Interaction between tomato and Fusarium oxysporum f. 
.•ip.lycopermi has been investigated in depth, becoming a model 
system for disease resistance response [9]. 
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Tomato mosaic virus (ToMV) is a positive-sense ssRNA virus 
belonging to the Tobamovirus genus. ToMV infects tomato plants 
systemically, causing mosaic symptoms, which are characterized 
by intermingled light and dark green regions [10]. The Tm-2 gene 
of tomato and its allelic form, Tm-2(2), confer resistance to 
Tomato mosaic virus (ToMV) and encode members of the CC- 
NB-ARC-LRR protein class [1 1]. 

These two pathogens have different infection strategies: ToMVis 
a biotrophic leaf pathogen that causes tomato leaf mold, whereas 
Fol is a saprotrophic soUborne vascular pathogen that challenges 
several Solanaccac". The aim of this study was to compare tomato 
global transcriptional profiles in response to host attack by ToMV 
and Fol in order to identify genomic differences and similarities in 
incompatible interactions between a foliar and a vascular 
pathogen, making use of the recently sequenced tomato genome. 
First we examined global tomato transcriptional profiling during 
both incompatible interactions to compare transcriptional changes 
occurring in tomato plants. Then we performed in-depth gene 
annotation, highlighting enriched GO categories and any meta- 
bolic perturbations arising. Finally, we explored the genome 
arrangement of expressed genes along the chromosome in order to 
connect genomic and transcriptional events. 

Materials and Methods 

Plant Material and Inoculation Protocol 

For our experiments we used the Rosso Delta tomato 
[S.ljcoper.skum) variety, resistant to Fol races 1 and 2 and ToAdV 
race 0. Tomato plants were inoculated with Fol strain ATCC 
16605 containing the Avrl and Avr2, (Plant Research Interna- 
tional, NL). Plantiets were infected at the stage of expanded 
cotyledons: roots were cut and dipped in a suspension at a 
concentration of 1 x 1 0** conidia/ ml, according to CPVO technical 
protocol TP/044/3 [12] (Experiment I). Inoculation with Tomato 
Mosaic Virus was carried out with the strain GM6s of ToMV (Plant 
Research International, NL), according to the method reported in 
tiie CPVO technical protocol TP/044/3 [12] (Experiment II). 
Plants with expanded cotyledons were inoculated with the sap 
obtained from infected desiccated tomato leaves. Viral transmis- 
sion was ensured by mechanical tissue abrasion caused by 
diatomaceous powder. 

Sample Collection 

Two days after treatment, infected and non-infected tomato leaf 
samples of two independent replicas of Experiment I and II were 
collected. Roots and leaves were removed from the plants, 
weighed and immediately frozen in liquid nitrogen and stored at 
— 80°C. RNA was isolated from whole plants using the Rneasy 
Plant Kit according to the manual instructions (Quiagen Valencia, 
USA). RNA sample concentration was determined using a 
Nanodrop photometer. An Agilent 2100 Bioanalyzer was used 
to check the quality of RNA. 

Chip Design and Microarray Hybridization 

Transcriptomic analysis was performed using a 90 K Toma- 
tArray 2.0 microarray synthesized using the CombiMatrix 
platform at the Plant Functional Genomics Center of the 
University of Verona. The chip (Tomato Array 2.0) carries 25,789 
non redundant probes (23,282 unique probes and 2,507 probes 
with more than one target) randomly distributed in triplicate 
across the array, each comprising a 35^0-mer oligonucleotide. 
The source of sequence information included tentative consensus 
sequences (TCs) derived from the DFCI Tomato Gene Index 
Release 12.0 and expressed sequence tags. Eight bacterial 



oligonucleotide sequences provided by CombiMatrix, eight probes 
designed on eight Ambion spikes and 40 probes based on Bacillus 
anthracis, Haemophilus ducreyi and Alteromonas phage sequences were 
used as negative controls. Complete description of the chip is 
available at the Gene Expression Omnibus under the platform 
accession GPL 13934. Microarray analysis was used to im-estigate 
tomato gene expression profiles two days after the infection with 
Fol and ToMV, comparing with the uninfected control profile. 
After checking the rjuantity and quality by spectrophotometry 
using NanoDrop 1 ()()() (Thermo Scientific), total RNA (2 |J.g) was 
amplified and labeled using the RNA am[)ULSe kit (Kreatech). 
For array hybridization, according to the manufacturer's recom- 
mendations, 4 |j,g of labeled aRNA was employed. Cy5 labeled 
aRNA was hybridized with the microarray at 45°C for 16 h; then 
arrays were washed with hybridization Solution at 45°C for 5 
minutes and 3xSS, 0.5xSSPET, PBST wash and PBS wash at 
room temperature. Pre-hybridization, hybridization, washing and 
imaging steps were carried out according to the manufacturer's 
instructions. Microarrays were stripped for reuse (up to three 
reuses per chip) with the CombiMatrix CustomArrayTM Strip- 
ping Kit according to the manufacturer's instructions (product 
number 610049). After hybridization and washing, the microarray 
was scanned using a Perkin Elmer Scan Array 4000XL (software 
ScanArray Express Microarray Analysis System v4.0). 

Data Analysis 

Gene expression levels corresponding to eight microarrays were 
processed using R software (R Core Team 2013 [13], URL 
http://www.R-project.org/.) and the limma package [14]. Raw 
data were investigated for quality assessment, preprocessed and 
normalized using a suitable subset of control probes available in 
the CombiMatrix array and a quantUe normalization technique. 
Three technical replicates within each array and two biological 
replicates were employed to assess differential expression for each 
experiment (Experiments 1 and 2) to compare the different 
experimental conditions (non-inoculated vs inoculated) using a 
linear model for microarray [15]. The significance of the 
differential expression was assessed taking into account the 
multiple testing setting and controlling the False Discovery Rate 
(FDR) at FDR = 0.01 and the empirical Bayes moderated t- 
statistics available in the limma package. Linear model fitting was 
used to assess differential expression of the two pairs of contrasted 
experimental conditions (non-inoculated vs inoculated), taking into 
account the presence of technical replicates as well as biological 
replicates. Differential expression data are deposited in Gene 
Expression Omnibus (GEO) database repository with accession 
number (GSE52336). 

Annotation of the Gene Chip Probes 

An in-house pipeline was developed to annotate tomato 
tentative consensus sequences (TCs) used to develop CombiMatrix 
CustomArrayTM probes. The queried tomato genes were 
identified by mapping TC sequences to the tomato CDS sequence 
using BLASTn (E-value le-3). The latest version of the tomato gfI3 
annotation files was parsed to extract the cds sequences of genes 
probed. Blast2GO http://blast2go.bioinfo.cipf es/was used to 
provide automatic high-throughput annotation, gene- ontology 
mapping and categorization of tomato proteins identified. An 
expectation value threshold of le-6 in BLASTp analysis was 
performed. Blast2GO was used for the statistical analysis of GO- 
term frequency differences. The enrichment analysis of the GO 
terms was based on Fisher's exact test and corrects for multiple 
testing. A cut-oflF FDR value of 0.05 was used. An interactive 
graph of the GO category enrichment REVIGO htt://revigo.org/ 



PLOS ONE I www.plosone.org 



2 



May 2014 I Volume 9 | Issue 5 | e94963 



Pathogen Infected Tomato Transcriptome Dynamics 



and a graphic representation Cytoscape htt://cytoscape.org// 
were also performed. 

Results 

Identification of Differentially Expressed Genes Induced 
by f. oxysporum f. sp. lycopersici (Fol) and ToMV 
Inoculation 

The transcriptional responses of resistant tomato seedlings, 
inoculated with Fol and ToMV, were evaluated by querying 
15,734 tomato genes. This analysis allowed us to compare changes 
occurring during both incompatible interactions with a vascular 
fungus and with a virus. DifiFerentiaUy expressed tomato genes 
(FDR-adjusted p-value <0.01) were identified at 2 DPI, compar- 
ing inoculated and non-inoculated plants. The time point was 
selected in order to identify genes involved in the initial stages of 
the defense against the pathogens. 

In Fol incompatible interaction, transcriptome variation resulted 
in 3,753 differentially induced genes. In particular, 2,392 genes 
(about 64%) were up-regulated (Fig. lA), indicating considerable 
gene activation during infection. As for the vascular pathogens, 
3,501 transcriptional changes were monitored in tomato upon 
inoculation with ToMV, of which about 2000 (52%) were 
overexpressed. When the total number of differentially regulated 
genes between the two incompatible interactions were compared, 
roughly half (2,205 genes) overlapped, of which a small subset of 
131 genes had the opposite expression direction during the two 
interactions. As shown in Fig. IB, a marked host gene induction 
was evident, since more than half of the differentially expressed 
genes appeared to be induced in botii interactions. Transcriptional 
changes during the two tomato-pathogen interactions were further 
investigated in order to find the best strategy for building a 
biological model for a plant-pathogen system. 

Gene Set Enrichment Analysis 

We performed a GO-term annotation analysis of all transcripts 
identified by gene chip probe matching. Through this analysis, we 
were able to assign functional annotations to 15,378 chip 
transcripts. In order to facilitate the formulation of biological 
hypotheses, a categorization of gene expression data was 
performed. Indeed, the Gene Ontology research effort was 
extremely useful for structuring data description. 

In the tomato- ToMF interaction 226 enriched GO terms were 
detected, 70% of which belonged to a biological process (P) 
category, while in the fungus interaction fewer enriched categories 
(185) were identified, 78% of which belonged to a biological 
process (P) category (Fig. 2A). Comparing the total number of GO 
categories between the two incompatible interactions, roughly 
60% overlapped. 

A heatmap obtained through hierarchical analysis evidenced 
four specific GO term clusters. Clusters 1 and 4 included GO 
ontology terms that did not show significant differences in two 
experiments. Cluster 2 contained about 90% of the enriched GO 
terms. Subcluster 2b consisted of specific GO terms of tomato-Fol 
interaction, whilst subclusters 2a and 2d as well as cluster 3 
included tomato- ToMF interaction specific ontology terms 
(Fig. 2B). 

To reduce redundancy of functional categories between the 
enriched GO terms, we used the semantic similarity approach. 
The gene ontology interactive graph-based network produced 
(Fig. 3) encapsulated functional homology between genes of the 
two tomato-pathogen interactions. In the tomato-Fol interaction, 
the network comprised 53 nodes and 310 edges, while for the 
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Figure 1. Comparison of tomato differentially regulated genes 
during interaction with Fusarium oxysporum f. sp. lycopersici (Fol) 
and Tomato iVIosaic Virus (ToMV) at 2 DPI post-inoculation. A, 

Venn diagram of differentially expressed genes in Fol and ToMV tomato 
incompatible interactions. B, Histogram of up regulated and down 
regulate genes in tomato-Fo/ interaction and tomato- FoMV interaction. 
doi:10.1371/journal.pone.0094963.g001 

tomato- ToMF interaction the network was more complex (58 
nodes and 383 edges). 

The network of tomato-ft/ interaction showed that the biotic 
stimulus node (GO:0009607) correlated with the five internal 
nodes with gene ontology terms related to external stimuli 
response (GO:0009637, GO:0009746, GO:0009628, 
GO:0042742, GO:0042744). The GO:0009607 category includes 
a large number of genes that play a role in signal transduction and 
regulation of gene expression in the defense response, representing 
about 10% of differentially expressed genes in tomato-Fo/ 
interaction. The external stimuli sub-network was also connected 
with hydrogen peroxide catabohsm node (GO:0042744). It 
contained 28 genes, 70% of which were overexpressed. In 
particular, four haem peroxidases (Solyc01g006290.2.1; So- 
Iyc01g006300.2.1; Solyc06g050440.2.1; Solycl lg018800.1.1), 
which play a role in host defense by inhibiting the hyphal 
extension of invading pathogens, were included in this node. 

Fig. 3B shows the network obtained with significant GO-term 
changes in tomato- ToAfF interaction. The biotic stimuli 
GO:0009607 node correlated with four GO-term nodes 
(GO:0009637, GO:0006950, GO:00 10038, GO:0009628), associ- 
ated to 264 genes. Interestingly, no external link was evidenced for 
this interaction. These data indicate that incompatible interactions 
induced both well documented stress-responsive genes as well as 
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Figure 2. Heatmap of Gene Ontology enriched terms in the two interactions. Gene Ontology terms enriched are obtained comparing a 
differential group (differentially regulated genes In each Interaction) whit a reference group (all gene of chip) using Fisher's Exact Test, with Multiple 
Testing Correction of FDR, estimated at P value of <0.05. A, The rows of heatmap represent GO-terms and the columns represent samples. Each cell 
Is colorized based on the number of genes associated with that category GO In that particular sample. B, Subclusters composed of specific GO-terms 
of the tomato-pathogen Interactions. The green boxes highlight specific GO terms of tomato-Fo/ Interaction, the red boxes specific GO terms of 
tomato-ToMl^ Interaction. 
dol:10.1371/journal.pone.0094963.g002 



unknown genes that might play a role in multi-stress responses. 
Investigation of single GO categories could evidence GO terms 
related to each assessed specific response. 

Investigation of Perturbed Biological Processes 

To obtain an overview of processes involved in fungus and virus 
incompatible reactions in-depth GO (gene ontology) term analysis 
was performed. Significant differences were observed in specific 
regulated biological processes between the two incompatible 
interactions. Notably, 30 specific enriched GO-term categories 
for tomaXo-Fol interaction, and 7 1 for tomato- ToA4V interaction 
were identified (Table SI and S2 in File SI). 

Tomato-Fo/ interaction revealed changes in cell structure, light 
perception, iron transport, aromatic amino acid synthesis and 
ROS response. In particular, we found six specific GO terms 



associated to homeostatic process (GO:0006873 cellular ion 
homeostasis; GO:0055080 cation homeostasis; GO:0050801 ion 
homeostasis; GO:0019725 cellular homeostasis; GO;0048878 
chemical homeostasis and GO:0042592 homeostatic process) that 
included 104 genes. This strongly suggests that homeostasis plays 
an important role in tomato defense to Fol. Interestingly, the 
master gene of inflammation, NF-kB (Solyc02g094530.1.1) was 
also up-regulated in tomato-i^o/ interaction. This gene is a key 
player in anti-apoptotic signaling and is able to prevent apoptotic 
signaling by inhibiting map-kinases. Indeed, tomato map-kinases 
(Solycl2g0 19460. 1.1; Solycl lg072630.1.1) were found to be 
heavily down-regulated during this interaction. Table 1 shows 
14 auxin-binding genes and two abscisic acid receptors up- 
regulated in tomato-Fo/ interaction. 
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Figure 3. The interactive graph-based enriched gene-ontology (A) in fungus-tomato interaction and (B) in virus-tomato interaction. 

doi:10.1371/journal.pone.0094963.g003 



Photosynthetic and carbohydrate derivative processes, nitrogen 
biosynthesis, therpene and carotenoid metabolism, and cadmium/ 
copper function were challenged in tomato- ToAfF interaction. 
Three specific GO terms associated to photosynthetic processes 
(GO:0009853 photorespiration; 00:0019685 photosynthesis dark 
reaction and GO:0009773 photosynthetic electron transport in 
photosystem I) were identified and were investigated firrther. 
Table 2 shows the enriched categories associated to photosynthetic 
processes containing 191, genes over 86% of which were down- 
regulated. Interestingly, photosynthetic gene transcription repres- 



sion was inversely correlated with the pathogenesis-related gene 
induction (Table 3). Moreover, a close linkage clearly emerged 
between the plant carbohydrate status and the resulting plant 
pathogen interaction. We found 1 7 carbon frxation enzymes 
differentially expressed during tomato- ToMF interaction. Further- 
more, several genes involved directly in the Calvin cycle were 
down-expressed (Table 4), and (3-fructofuranosidase (So- 
Iycl0g085640. 1. 1), a key invertase involved in sink source 
portioning, was overexpressed. We also found up-regulated six 
gibberellin-modulated genes and two genes involved in SA 



Table 1. List of abscisic acid receptors, auxin response factors (ARFs) and auxin response genes up-regulated in the tomato-Fo/ 
interaction. 



Gene ID 


Probe ID 


Putative function 


GO-ID terms 


Solyc01g091 030.2.1 


TCI 96630 


Auxin responsive SAUR protein 


GO:0009733; GO:0009862; GO:0003674 


Solyc02g077560.2.1 


TCI 92465 


Auxin response factor 


GO:0009725: GO:0009734: GO:0009850: GO:0009733 


Solyc03g007310.2.1 


TCI 97392 


Abscisic acid receptor 


GO:0010427 


Solyc03g031 970.2.1 


TC207784 


Auxin response factor 


GO:0009734 


Solyc03g120500.2.1 


TCI 99757 


Auxin responsive protein 


GO:0009734 


Solyc04g074980.2.1 


TC207543 


Auxin F-box protein 


GO:0009734 


Solyc05g0474e0.2.1 


TC205766 


Auxin response factor 


GO:0009734 


SolycOegOeil 80.1.1 


TCI 93498 


Abscisic acid receptor 


GO:0010427: GO:0009738 


Solyc07g01 61 80.2.1 


TC207832 


Auxin response factor 


GO:0009734 


Solyc07g043ei 0.2.1 


TC207520 


Auxin response factor 


GO:0003677: GO:0009734 


Solyc07g0e3850.2.1 


TC2 16744 


GH3 auxin-responsive promoter 


GO:0010252: GO:0009734 


Solyc09g01 4380.2.1 


TC2 16697 


Auxin transporter-lil<e protein 


GO:0009734: GO:0010328: GO:0060919: GO:C009926: GO:0010011 


Solyc10g050710.1.1 


TC202297 


GH3 auxin-responsive promoter 


GO:0010279: GO:0009733 


Solyc10g083320.1.1 


TCI 94758 


Auxin responsive SAUR protein 


GO:0009733 


Solyc11g0691 90.1.1 


TC208143 


Auxin response factor 


GO:0009734 


Solyc12g005310.1.1 


TC203007 


GH3 auxin-responsive promoter 


GO:0010279: GO:0010279: GO:0009734 


For each entry is reported: Gene and probe 
doi:l 0.1 371/journal.pone.0094963.t001 


ID, description of putative function and GO- 


ID term correlated. 
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Table 2. List of the GO-terms involved in the photosynthetic process during the tomato-ro/MI/ interaction. 



GO ID 


Term 


Category 


n. genes 


% downeregulated genes 


1 


photosystGm 


cellular component 


25 


1 00 




photosystem 1 


cellular component 


1 6 


1 00 




photosystem II 


cellular component 


1 8 


1 00 


ijl_I.UUU7/03 


photosynthesis, light harvesting 


biological process 


1 7 


1 00 


GO:0009767 


photosynthetic electron transport chain 


biological process 


24 


100 


GO:0009773 


photosynthetic electron transport in photosystem 1 


biological process 


15 


86,7 


GO:00 10207 


photosystem II assembly 


biological process 


50 


94 


GO:0015979 


photosynthesis 


biological process 


125 


91,2 


GO:00 19684 


photosynthesis, light reaction 


biological process 


93 


94,6 


GO:00 19685 


photosynthesis, dark reaction 


biological process 


12 


91,7 


GO:0034357 


photosynthetic membrane 


cellular component 


127 


89,8 



For each GO- ID term, the category, the number of genes identified and the percentage of down regulation revealed is reported. 
doi:l 0.1 371 /journal.pone.0094963.t002 



(Solyc01g014320.2.1) and JA (Solyc07g042 170.2.1) synthesis 
(Table 5). 

Relationship between Genomics and Transcriptional 
Changes 

In order to investigate the chromosomal arrangement of 
differential expressed sequences identified by tomato chip, a 
genomic expression map was constructed (Fig. 4). Our results 
indicated that differentially expressed genes (about 33% of all 
genes analyzed) during the two interactions were arranged along 
the chromosomes with different levels of expression. The fungus 
and virus distribution showed slight differences. Indeed, the 
number of expressed genes identified for single chromosomes was 
always greater in tomato- ToMF interaction than in tomato-ft/ 
interaction, with the exception of chromosome 8. Chromosomal 
areas (on chromosomes 1, 2, 8, 11 and 12) showed a higher density 
of up-regulated genes in the fungus-interaction. By contrast, on 
chromosome 8 we found a region with a higher density of down- 
regulated genes in the virus-interaction. In both interactions, the 
genomic regions rich in differentially expressed genes were located 
in the vicinity of telomeres, inter alia the richest regions of genes 
were the long arm of chromosomes 10 and 12. 

The genome-wide distribution of differentially expressed genes, 
based both on chromosome size and gene density was not 
uniformly distributed along the genome (chromosome size: fungus- 
interaction x2 = 107; virus-interaction x2 = 117. Gene density: x2 
value included between 27 and 170). Fig. 5 and 6 showed the 
genomic distribution of under-expressed and over-expressed genes 
during the two interactions. If we use the genome location of the 
coding genes spotted on the microarray as a reference we could 
highlight a significant difference in the distribution of over- 
expressed genes as well as under-expressed ones (Kolmogorov- 
Smirnov and Cramer von Mises tests in Table S3 and S4 in File 
SI). Indeed, differentially expressed genes (p-value <0.10) were 
evidenced on all chromosomes, except for chromosomes 2, 4 and 
7. Stronger evidence of peculiar chromosomal displacement of 
over-expressed genes may be found on chromosomes 3, 5, 6 and 1 
(p-value <0.05). Significantly differential displacement of under- 
expressed genes (p-value <0.05) was detected mainly in chromo- 
somes 5 and 9. This suggested the presence of a dynamic genomic 
landscapes involved in pathogen response. 



For each interaction studied, distribution of expressed specific 
genes is shown in Fig. 7. Most were located in the same regions, 
but clear differences could be highlighted. The 151 common 
differentially expressed genes with opposite direction were mainly 
arranged on chromosomes 7, 9 11 and 12. On chromosomes 9 
and 11, in the vicinity of genes Tm2(2) conferring resistance to 
Tomato mosaic virus and I and 12 conferring resistance to Fusarium 
oxy.spomm, we highlighted differences in gene expression between 
the two experiments. Eleven genes were specifically expressed in 
ToMV interaction in this area of chromosome 9 (including: 
pentatricopeptide repeat, ATP dependent RNA hehcase, cyto- 
chrome P450 and pectinesterase) and 19 of the 12 region of 
chromosome 1 1 in Fol interaction (including: BEL 1 -like homeo- 
domain protein 8, Calmodulin-lilie protein and MYB transcription 
factor and two proteins of unknown function).whereas around of I 
locus in the 1 Mb region delimitated from TG523 and CT168 
[16], 8 specific genes were identified (including a Pumilio-like 
IPRO 11989 protein and Serine-threonine protein kinase). Fig. 8 
shows differential expressed pathogen recognition genes (NB- 
ARC-LRR, RLP, RLK) both in tomato-ft/ and tomato- ToM7 
interactions. Pronounced activation of all pathogen recognition 
gene classes assessed was evidenced in tomato-Fo/ interaction since 
70% of all the differentially expressed genes were up-regulated. 
Instead, in the tomato- ToAfF interaction the number of up- 
regulated and down-regulated pathogen recognition genes was the 
same: 45% of all the differentially regulated pathogen recognition 
genes overlapped. Of these genes, 10 had the opposite expression 
in the two experiments (Table S5 and S6 in File S2). 

Discussion 

Plant-expressed molecules during pathogen challenging can give 
insights into the underlying defense mechanisms. Plants have 
evolved a defense system against microbial pathogens that involves 
the regulation of gene expression, cascade signaling activation, 
hormone balancing and synthesis of defensive metabolites [17]. 
Several microarray experiments have been successfully carried out 
in the sphere of molecular plant-microbe interaction, shedding 
light on the mechanisms controlling plant disease resistance and 
the crosstalk among the signaling pathways involved [4]. More 
than 3500 genes, out of 15,734 assessed, were activated in tomato 
ft/ and ToAfF interactions, more than 60% of which overlapped. 
Pathogens from different kingdoms deploy independently evolved 
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Table 3. List of up-regulated pathogenesis-related proteins involved in tomato-IoMI/ interaction. 



Protein ID 


Probe ID 


PRs family 


Enzyme 


^Fold change 


SolycOl g006290.2.1 


TCI 97719 


PR-9 


Peroxidase 


5.5 


SolycOl g067860.2.1 


TCI 95695 


PR-9 


Peroxidase 


2.3 


SolycOl g 108320.2.1 


TC209477 


PR-9 


Peroxidase 


2.3 


Solyc02g064970.2.1 


TC201774 


PR-9 


Peroxidase 


8.9 


Solyc02g077300.1 .1 


TC20541 1 


PR-9 


Peroxidase 


1.6 


Solyc02g0795 10.2.1 


TCI 96561 


PR-9 


Peroxidase 


5.7 


Solyc02g080530.2.1 


TC209879 


PR-9 


Peroxidase 


6.9 


Solyc02g082960.2.1 


TC200207 


PR-3 


Endochitinase 


1 .7 


Solyc02g084780.2.1 


TC201 521 


PR-9 


Peroxidase 


6.3 


Solyc02g0941 80.2.1 


TC21 1 238 


PR-9 


Peroxidase 


2.6 


Solyc03g005200.2.1 


TC203320 


PR-14 


Non-specific lipid-transfer protein 


2.8 


Solyc03g03371 0.2.1 


TC203104 


PR-9 


Peroxidase 


7.1 


Solyc03g0441 00.2.1 


TCI 99639 


PR-9 


Peroxidase 


5.8 


Solyc04g007750.2.1 


TCI 99935 


PR-10 


Major latex-like protein 


3.2 


Solyc04g071 890.2.1 


TCI 93 192 


PR-9 


Peroxidase 


4.3 


Solyc04g072000.2.1 


TC2 12476 


PR-3 


Chitinase 


1.2 


Solyc04g080760.2.1 


TC206568 


PR-9 


Peroxidase 


2.1 


Solyc05g04601 0.2.1 


TC202014 


PR-9 


Peroxidase 


2.4 


Solyc05g0501 30.2.1 


TCI 94980 


PR-1 1 


Acidic chitinase 


2.8 


Solyc06g050440.2.1 


TC21 1424 


PR-9 


Peroxidase 


6.9 


Solyc06g072220.1 .1 


TC2 10670 


PR-6 


Kunitz trypsin inhibitor 


SS 


Solyc06g076630.2.1 


TC201658 


PR-9 


Peroxidase 


7.9 


Solyc07g005380.2.1 


TCI 94633 


PR-10 


Norcoclaurine synthase 


2.6 


Solyc07g009260.2.1 


TC2 10032 


PR-1 2 


Defensin-like protein 


2.3 


Solyc07g01 7880.2.1 


TC20971 0 


PR-9 


Peroxidase 


3.8 


Solyc07g0525 10.2.1 


TC2081 95 


PR-9 


Peroxidase 


8.7 


Solyc08g067510.1.1 


TC213655 


PR-14 


Non-specific lipid-transfer protein 


1.5 


Solyc09g007520.2.1 


TC200697 


PR-9 


Peroxidase 


3.9 


Solyc09g065430.2.1 


TC2 14938 


PR-14 


Non-specific lipid-transfer protein 


2.1 


Solyc09g082280.2.1 


TC209373 


PR-14 


Non-specific lipid-transfer protein 


8.9 


Solyc10g076220.1.1 


TC209457 


PR-9 


Peroxidase 


6.4 


Solyc11g072760.1.1 


TC215379 


PR-11 


Chitinase a 


3.7 


Solyc12g01 7660.1.1 


TCI 98426 


PR-9 


Peroxidase 


2.9 



^Fold change is calculated using the signal log2 ratio. 

For each entry it is reported: gene and probe ID, the PRs-family to which belongs and the description of its function. 
doi:1 0.1 371 /journal.pone.0094963.t003 



virulence protein.s that interact with a hmited set of highly 
connected plant cellular hubs [18]. Network analysis suggests that 
plant response to Fol and ToMV infection involved activation or 
repression of genes implicated in pathways shared by most 
responses to environmental stimuli. However, GO enrichment 
analysis evidenced specific enriched categories in both interac- 
tions. Response to ToMV seemed more multifaceted, since more 
than 70 specific categories were enriched versus the 30 detected in 
Fol interaction. Biotrophic plant pathogens are generally accepted 
to have a more intricate biological interaction with their host plant 
than saprophytic plant pathogenic fungi [19]. Once an attack is 
perceived, plant metabolism must balance potentially competing 
demands for resources to support defense versus requirements for 
ceUular maintenance, growth and reproduction [20;21;22;23]. 

In tomato-Fo/ interaction our investigation evidenced a number 
of overexpressed genes associated to maintenance of cellular 



structures and cellular homeostasis. These are very important 
metabolic activities required by plants to survive fungus-inflicted 
stresses. Fusarium oxjsporum is a necrotrophic fungus [24] that kills 
host cells prior to the infection through the predicted deployment 
of toxins and enzymes, which induces cell death. Hence the 
expression of anti-apoptosis genes could confer resistance to Fol 
[25]. The NF-kB signaling enhancement found in this work can 
inhibit apoptosis and combat the effects of oxidative stress by 
switching down the anti-apoptotic signals mediated by the map- 
kinase cascade. These events lead us to hypothesize that the 
resistant plant challenged by Fol induces a homeostatic response to 
prevent the fungus from killing the cells to obtain its nourishment. 
Fourteen auxin-binding genes were also activated in this 
interaction. Auxin signaling is an important phytoregulator for 
resistance to saprophytic/necrotrophic fungi [26]. The degrada- 
tion of AUX/IAA proteins allows activation of auxin response 
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Table 4. List of genes involved in carbon fixation differentially expressed during the tomato-Io/MI/ interaction. 



Gene ID 


Function Description 


Ezyme ID 


Enzyme class 


Up/Down-expressed 


SolycOSgOl 3860.2.1 


Malic enzyme 


ec:1.1.1.39 


dehydrogenase 


Up 


Solyc12g044600.2.1 


Malic enzyme 


eel. 1.1. 40 


dehydrogenase 


Down 


Solyc03g071 590.2.1 


Malate dehydrogenase 


ec:1.1.1.82 


dehydrogenase 


Up 


Solyc04g009030.2.1 


Glyceraldehyde-3-phosphate dehydrogenase 


eel. 2.1. 13 


dehydrogenase 


Down 


Solyc12g094640.1.1 


Glyceraldehyde-3-phosphate dehydrogenase 


ec:1.2.1.13 


dehydrogenase 


Down 


Solyc10g018300.1.1 


Transketolase 


ec:2.2.1.1 


glycolaldehydetransferase 


Down 


Solyc05g050970.2.1 


Transketolase 


ec:2.2.1.1 


glycolaldehydetransferase 


Up 


Solyc08g079750.2.1 


1 -aminocyclopropane-1 -carboxylate synthase 
Aminotransferase 


ec:2.6.1.1 


transaminase 


Down 


Solyc07g032740.2.1 


Aspartate aminotransferase 


ec:2.6.1.1 


transaminase 


Ud 


Solyc07g05521 0.2.1 


Aspartate aminotransferase 


ec:2.6.1.1 


transaminase 


Ud 


Solyc11g044840.1.1 


LL-diaminopimelate aminotransferase 


ec:2.6.1.1 


transaminase 


Ud 


SolycOSgOl 3380.2.1 


Alanine aminotransferase 


ec:2.6.1.2 


transaminase 


Down 


Solyc06g063090.2.1 


Alanine aminotransferase 


ec:2.6.1.2 


transaminase 


Ud 


Solyc08g076220.2.1 


Phosphoribulokinase 


ec:2.7.1.19 


phosphopentokinase 


Down 


Solyc04g008740.2.1 


Pyruvate kinase 


ec:2.7.1.40 




Ud 


Solyc07g066600.2.1 


Phosphoglycerate kinase 


ec:2.7.2.3 


kinase 


Down 


Solyc07g06661 0.2.1 


Phosphoglycerate kinase 


ec:2.7.2.3 


kinase 


Down 


Solyc09g011 81 0.2.1 


Fructose-1 6-bisphosphatase 


ec:3. 1.3.1 1 


hexose diphosphatase 


Down 


Solycl 0g086730.1 .1 


Fructose- 1 6-bisphosphatase 


ec:3. 1.3.1 1 


hexose diphosphatase 


Down 


Solyc12g0142S0.1.1 


Phosphoenolpyruvate carboxylase 


ec:4.1.1.31 


carboxylase 


Down 


Solycl Og007290.2.1 


Phosphoenolpyruvate carboxylase 


ec:4.1.1.31 


carboxylase 


Ud 


SolycOl g007330.2.1 


Ribulose bisphosphate carboxylase 


ec:4.1.1.39 


carboxylase 


Down 


Solyc03g034220.2.1 


Ribulose bisphosphate carboxylase 


ec:4.1.1.39 


carboxylase 


Down 


Solyc06g061 280.2.1 


Cinnamoyl-CoA reductase- like protein 


ec:4.1.1.49 


carboxykinase 


Up 


Solyc01g110360.2.1 


Fructose-bisphosphate aldolase 


ec:4.1 .2.1 3 


aldolase 


Down 


Solyc07g065900.2.1 


Fructose-bisphosphate aldolase 


ec:4.1.2.13 


aldolase 


Down 


Solyc09g009260.2.1 


Fructose-bisphosphate aldolase 


ec:4.1.2.13 


aldolase 


Up 


Solycl Og083570.1.1 


Fructose-bisphosphate aldolase 


ec:4.1.2.13 


aldolase 


Up 


Solyc03g1 15820.2.1 


Ribulose-5-phosphate-3-epimerase 


ec:5.1.3.1 


epimerase 


Down 


Solyc06g005490.2.1 


Triosephosphate 


ec:5.3.1.1 


isomerase 


Down 


For each gene is reported: the functional description, the enzyme ID, the enzyme class 
doi:1 0.1 371 /journal.pone.0094963.t004 


and direction of 


expression regulation. 





factors (ARFs) and the expression of auxin-responsive genes [27]. the sites of pathogen ingress constitutes one of the main host 

These protein families can liave an enormous number of factors that determine plant resistance to Fusarium wilt [29]. 

interactions, capable of fine-tuning specific responses within the Essential regulators of Jasmonate-ZIM domain proteins (JAZs) in 

auxin signaling pathway [28]. Moreover, the buildup of lAA near tomato-ft/ interaction were activated, suggesting attenuated JA 



Table 5. List of the gibberellin modulated genes, up-regulated 


in tomato-roMV interaction. 




Gene ID 


Probe ID 


Putative function 


GO-ID term 


Solyc02g083880.2.1 


TC203352 


Gibberellin-regulated protein 


GO:0005515 


Solyc0eg008870.2.1 


TC202708 


Gibberellin receptor 


GO:0005515; GO:0008152 


Solyc06g035530.2.1 


TC202928 


Gibberellin 20-oxidase-2 


GO:0016491; GO:0045544; GO:0055114 


Solyc10g005360.2.1 


TC204547 


Gibberellin 2-beta-dioxygenase 


GO:0016491,GO:0045543 


Solyc10g050880.1.1 


TCI 98803 


Gibberellin receptor 


GO:0004091 ,GO:00081 52 


Solyc11g072310.1.1 


TC211262 


Gibberellin 20-oxidase-3 


GO:001 6491 ,GO:0045544,GO:00551 1 4 


For each gene is reported: genes and probe 
doi:l 0.1 371/journal.pone.0094963.t005 


ID, the description of putative function and the GO-ID terms correlated. 
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Figure 4. Physical genome locations of tomato (Solanum lycopersicum) differentially regulated genes. A, Phiysical map of the tomato-Fo/ 
interaction. B, Physical map of the tomato-7o/Ml/ interaction. Tomato chromosomes are represented as blacl< horizontal bars, the approximate 
location of each gene is designated with vertical lines on chromosomes. The color used for each gene indicates specificity (red: differentially 
regulated genes in tomato-fungus interaction, blue: differentially regulated genes in tomato-virus interaction, green: overlapping differentially 
regulated genes. The cyan boxes on chromosomes 9 and 1 1 indicates the chromosomal regions in the vicinity of the Tm2(2) and 12 genes. The violet 
boxes on chromosomes 8, 10 and 12 indicates the richest regions of genes. 
doi:1 0.1 371 /journal.pone.0094963.g004 



signaling to reduce the senescence process. Indeed, also a pair of 
abscisic acid receptors were up-regulated in Fol interaction. 
Interactions between multiple components of ABA and the JA- 
ethylene signaling pathways modulate defense and stress-respon- 
sive gene expression in response to Fol, confirming that these 



hormones mediate resistance in necrothophic interaction [30;31]. 
The SA pathway proved to be switched down. Indeed, silencing a 
tomato gene encoding SA methyltransferase was found to enhance 
resistance to F. oxysporum i. sp. lycopmici biotic and abiotic stresses 
[32]. 
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Figure 8. Comparison of differentially regulated tomato pathogen recognition genes during interaction with Fusarium oxysporum 

(Fol) and Tomato Mosaic Virus (ToMV). Histograms display the differentially regulated pathogen recognition genes for each structural class (NL: NB- 
ARC-LRR proteins RLP Like Proteins; RLK Receptor-Like Kinases; Unicnown: genes that encode novel domain associations or single domains) during 
the two incompatible interactions. 
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In tomato- ToA/F interaction, we detected a reduced photosyn- 
tlietic activity and elevated carbohydrate catabolism. Pathogen 
attack may liave caused metabolic alterations that induced a 
hormone response for reprogramming the cellular metabolism. 
Tlie photosynthesis and assimilatory metabolism were switched off 
to initiate respiration and other processes required for defense. 
This occurrence reflected the reallocation of plant metabolites 
from normal growth processes to defensive functions after the 
elicitation of induced plant responses by virus infection [33]. A 
close linkage between the plant carbohydrate-status with the 
resulting plant pathogen interaction was also evident. During 
interaction with the virus, the host produced salicyhc acid and 
increased the production of an invertase enzyme that is able to 
redirect the flux of carbohydrate acquisition. The infection of 
ToMV could have led to the alteration of the concentration 
gradient of sucrose in phloem. The plant overexpressed P- 
fructofuranosidase (Solycl0g085640.1.1) to ensure the sucrose 
transport of source ceUs to the sink ceUs. Overexpression of 
invertase in tomato-virus interaction could be elicited by the 
expression of pathogenesis-related proteins and by salicylic acid in 
order to increase resistance to virus infection [34]. In order to 
ensure the success of defense, this mechanism appears to be 
crucial. The phenomenon of "high sugar resistance" was 
described long ago [35] and the finding that various pathogen- 
esis-related genes were sugar-inducible [36] supports this hypoth- 
esis. We also found six gibbereUin, one SAM dependent carboxyl 
methyltransferase (SolycOlgO 14320.2.1) and one jasmonate-zim 
domain (Solyc07g042 170.2.1) genes up-regulated, suggesting that 
GAs activated tomato immune responses to ToMV by modulating 
the levels of salicylic acid and/or jasmonic acid [37]. 



Plants rely on the innate immunity and on systemic signals 
emanating from metabolic alterations [38] . The ETI system seems 
to prompt response in the right direction thanks to metabolic clues 
and hormone signaling. Salicylic acid (SA) primarily triggers 
resistance against biotrophic pathogens, whereas a combination of 
jasmonic acid (J A) and ethylene (ET) signaling activates resistance 
to necrotrophic pathogens [39;40]. The metabohc changes 
associated with defense response in two incompatible tomato- 
pathogen interactions suggested that the response to the specific 
metabolic alteration (photosynthesis/carbohydrate metabolism 
and homeostasis activity) in tomato was pathogen-specific and 
contributed substantially to monogenetic gene resistance. Our 
results confirm that resistance to pathogens depends on a 
sophisticated interplay among different biological pathways and 
that hormonal directionality is critical to the outcome of a response 
(carbohydrate biosynthesis and homeostatic activity). Several 
hormones involved in pathogen perception, activation of defense 
products restricting pathogen invasion, have been identified [2]. 
How they trigger ETI defense is currendy unclear. In the absence 
of a pathogen, NB-ARC genes are in an autoinhibition state, 
which is relieved upon pathogen perception [41]. Although 
difierent forms of plant immunity share the same signaling 
mechanisms, they use the same mechanisms in very different ways 
[42]. The mechanisms that lead to rapid metabohsm switch and 
connection among all the defense pathways is still not clear. 
Signaling able to fine-tune the defense mechanism could be 
activated by genes in proximity of the main pathogen receptors. 

Genome organization of functional gene networks to tolerate 
alterations can result in plasticity. In two specific interactions the 
chromosome distribution of expressed genes showed wide 
overlapping regions, except for the region holding genes 12 and 
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Tm2(2). Indeed, genome regions can be enriched in genes with the 
specific fimction for fine-tuning gene expression in a compensatory 
way. In Fol interaction on chromosome 1 1 region harboring 12 gene 
a calmodulin gene, a myb factor and a BEL-1 like homeodomain 
protein 8 were specifically expressed. Calmodulin (CaM) plays an 
important role in sensing and transducing changes in cellular Ca2^ 
concentration in response to several biotic and abiotic stresses, In 
ToMV interaction Pentatricopeptide-repeat, ATP-dependent and 
RNA hehcase could be involved in the process of silencing and 
virus replication. 

A high level of NB-ARC genes activation and of other gene 
classes potentially involved in pathogen recognition was also 
found, suggesting a host-coordinated reaction of defense machin- 
ery to monitor integrity of cellular proteins. The NB-ARC system 
minimizes the cost of defense for the plant, as multiple NB-ARC- 
LRR proteins can be maintained at a low level in the absence of a 
pathogen and rapidly induced under pathogen attack through 
miRNA regulation [43] . Pathogen-encoded suppressors of RNA- 
silencing mechanisms might result in the induction of multiple NB- 
ARC-LRR defense proteins [44]. Investigation of dilferentially 
regulated pathogen recognition genes could lead to the identifi- 
cation of specific modulation patterns. 

Albeit still fragmented, our depirtion pro\ idcs a global view of 
the I-I2 and Tm2(2) mediated resistance process, in which the gene 
expression network may be considered the starting point to 
construct a genomic model for R-mediated response in the two 
pathosystems investigated. 
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